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Abstract 

We evaluate the amplitude for 77 — > 7r + 7r~ to two loops in the framework of 
chiral perturbation theory. The three new coupling constants that enter the 
result at this order in the low-energy expansion are estimated via resonance 
saturation. We discuss in addition the crossed channel processes 77r ± —* 
- in particular the charged pion polarizabilities - to the same accuracy. The 
predictions are compared with available experimental information. 
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1 Introduction 



Pion-pair production by photons at low energies may be used to test chiral perturbation 
theory (CHPT) CJ-pI. For neutral pions, the chiral calculation at leading order - generated 
by one-loop diagrams |], - is not in good agreement with the available data from 
Crystal Ball H even near the threshold. Recently, the next-to-leading order contribution, 
generated by two-loop diagrams, has been determined in Ref. ||. The corresponding 
enhancement of the cross section brings theory and experiment into agreement within the 
uncertainties of the calculation and the available data |J . 

In the process 77 — ► 7r + 7r~, the Born term (that is absent in the neutral channel) 
dominates the behavior of the cross section at low energies. Next-to-leading corrections 
(one-loop graphs) change the result very little ||, and the comparison with the available 
data from Mark II pi works quite well. The dispersive analysis performed in [11, 12] 



indicates how additional (multiloop) effects required by unitarity may influence the am- 
plitude. The authors find only a slight change as compared to the one-loop result for 
center-of-mass energies y/s < 700 MeV. It is one of the purposes of the present article to 
investigate the effect of the two-loop contributions to the amplitude. 

The Compton amplitude jir^ — > 771"^ may be obtained from 77 — > ir + ir~ by crossing. 
The leading term in the Taylor expansion of the Compton amplitude at threshold in powers 
of photon energies measures the charge of the pions, whereas the next-to-leading terms are 
parametrized through the electric (a n ) and magnetic (f3 n ) polarizabilities. These probe the 
rigidity of the charged pion against an external electromagnetic perturbation. At leading 



order in the chiral expansion, the result is [14-Q a n = —(3^ = 2.7 ± 0.4 q, With the 



two-loop Compton amplitude at hand, it is then straightforward to evaluate a w and p w at 
next-to-leading order in the quark mass expansion. The result of this analysis is presented 
below. Turning to the experimental aspect, we note that the experimental determination of 



the charged pion polarizabilities is not an easy task. At Lebedev [|T7fl and at Serpukhov [18] , 
relevant experiments have been performed some time ago. Unfortunately the results of 
these measurements do not agree very well among each other, such that a comparison with 
the theoretical prediction is not conclusive. In order to clarify the situation, additional 
experiments are planned [|i~9|]. 

The article is organized as follows: After setting up the kinematics and illuminating 
some aspects of CHPT, the effective lagrangians £2,^4 and Cq are discussed. For com- 
pleteness we display in section [| the tree- and one-loop result worked out earlier in ||. 
We then evaluate all diagrams occurring at 0(p & ), including also the pion mass and the 
wave function renormalization constant at two-loop order (sections |(| 0) . The amplitude 
for 77 — > 7r + 7r~ is displayed in section B that also contains a numerical discussion of the 
corresponding cross section. Section || is devoted a discussion of Compton scattering and 
pion electric and magnetic polarizabilities at two-loop order. Finally, a summary and 
concluding remarks are given in section |l0| The notation used in the text is collected in 
appendix |A[ 



2 We express the polarizabilities in units of 10 4 fm 3 throughout. 
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2 Kinematics 

The amplitude describing the process 77 — > 7t + 7t~ may be extracted from the matrix 
element 

( ir + ( Pl ),ir-( P2 ) out I 7fe, Ai), 7(92, A 2 ) in ) = i(2vr) 4 5 4 (P / - P) T^ +w ~, (2.1) 
where 

r77 _- +7r - = e 2 e^(g 1 ,A 1 )6^( (72 ,A 2 )W;r^^ , 

W^T***" = » / dxe"*^*^ ( 7r+( Pl )7r-(p 2 ) out \ T j„(x)j u (y) | ) . (2.2) 

denotes the electromagnetic current and a = e 2 /4-7T ~ 1/137 is the electromagnetic cou- 
pling. It is convenient to change the pion coordinates according to (71^,71°) — * (tt 1 , tt 2 , 7r 3 ) 
and instead of 7r + 7r~-production, we consider in the following the process 77 —* tt 1 !! 1 , 
with 

W ry— + -- = -Wrf-™ 1 * 1 = -V° , (2.3) 

where the relative minus sign stems from the Condon-Shortly phase convention. The 
correlator V^ u may be decomposed into four Lorentz invariant amplitudes A c ,B C , C° ,D° 
(see e.g. |J), 

V° v = A c {s,t,u)T 1 ^ + B c (s,t,u)T 2 ^ + C c {s,t,u)T^ v + D c {s,t,u)T^ u , 

rp 8 

1 lfj,u — ifiixv — 1li/l2fj, i 

T 2liv = 2sA^A u - v 2 g liv - 2v {q\ u A^ - q 2f ,A u ) , 
Ts^v = qifj,q2v 7 

T^v = s (qi^A u - q 2v A^) - v (qx^qiv + q2^2v) , 
A„ = (Pi-P2) M) (2.4) 

where 

s = (qi +q2? =4q 2 , 

t = (pi-qif = M 2 -2q 2 (I- (3(s) cos 6) , 
u = ( P2 -qi) 2 = M 2 -2q 2 (l+P(s)cose) , 

v = t-u, (2.5) 

are the standard Mandelstam variables, f3(s) = (1 — 4M 2 /s) 1 / 2 is the velocity of the 
produced pions and 6 denotes the scattering angle, q\.p\ = |<fi||pi| cos# (all quantities in 
the center-of-mass system). The tensor V^ v satisfies furthermore the Ward identities 

<£V%, = fiV° = . (2.6) 

The amplitudes A° and B c are analytic functions of the variables s, t and u, symmetric 
under crossing (t,u) — ► (u,t). The remaining quantities C c and D° do not contribute to 
the cross section (gauge invariance). 
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The differential cross section for unpolarized photons in the center-of-mass system 
reads with the normalization ( pi | p2 ) = 2(2ir) 3 p ( [5 3 (pi — p^), 

H c (s,t) = |ff^ + | 2 + |^_| 2 . (2.7) 

Here we have introduced the helicity amplitudes H+± corresponding to photon helicity 
differences A = 0, 2, respectively. They are defined in terms of the amplitudes AF and B c , 

Hl + = A c + 2(4M 2 - S )B C , 

S 

and have partial wave expansion involving even J > A, 

H ++ = E hl + ( S )di(cos8) , 

,7=0,2,4... 

H%- = E /»£_M4(«>s0). (2.9) 

J=2,4,6... 

In the text we also use the symbols 

frC _l_ ttC _ ttC 
n ++ — IJ ++ n B++ > 

H%_ = H c + _-H c B+ _, (2.10) 

where -ff^+i denote the helicity amplitudes generated by the Born contributions. 

For comparison with experimental results, it is convenient to present also the total cross 
section for the case where | cos 6 \ is less than some fixed value Z , 

a(s, | cos 6 | < Z) = ^ f + dt H c (s, t) (2.11) 

with 

t± = M*-2q 2 (lTP(s)Z) . (2.12) 



3 Low energy expansion 

We consider QCD with two flavours u and d in the isospin limit m u = = rh. Since the 
masses of these two quarks are small, the deviation from chiral symmetry may be studied 
by treating the quark mass term as a perturbation using the framework of an effective chiral 
lagrangian (CHPT) |l]-|5|. This effective lagrangian relies only on the solid assumptions 
of spontaneously broken SU(2)l X SU (2)r chiral symmetry, Lorentz invariance and low 
momentum transfer. It can be expanded in a series of terms with increasing numbers of 
derivatives and powers of quark masses, 

C eff = C 2 + U + C & + ..., (3.1) 

where C^n) denotes a term of order p 2n . Each of these terms contain a couple of monomials 
with couplings which are not fixed by symmetry constraints, but they are fixed by the 



4 



dynamics of the underlying theory through the renormalization group invariant scale Aqcd 
and the heavy quark masses. Up to now it has not been possible to reliably calculate 
them directly from the QCD lagrangian. The best way to determine their values is by 
comparison with experimental data. Green functions of the theory are generated by the 
vacuum-to-vacuum amplitude pLB 



iZ(v,a,s, P ) = { | Q . n ) = I [du]el Jd^C eff (U,v,a,s,p) ^ ^ 



6 — \ ^out I u in / v,a,s,p 

where the dynamical variables (pions) are contained in the matrix U. The generating 
functional Z admits an expansion in powers of the external momenta and the quark 
masses 

Z(v, a,s,p) = Z 2 + Z 4 + Z 6 + . . . , (3.3) 

where Z n denotes a term of order p n . The contribution of loops to the generating functional 
is suppressed with respect to tree diagrams. At order p 2 the contribution to the amplitude 
stems from £2 only via tree diagrams. At next -to-leading order one has to evaluate one- 
loop graphs generated by £2 and adding the tree graphs from £4. Diagrams at C(p 6 ) 
contribute via Z§ from £2 + £4 + Cq with tree, from £2 + £4 with one-loop and from £2 
with two-loop diagrams. The one-loop diagrams at 0(p 6 ) contain one vertex from £4. 
We write the corresponding expansion of the amplitude as 

I = I 2 + h + I 6 + ..., I = V^,A C ,B C , (3.4) 

where it is understood that Z n generates I n . For comprehensive reviews on chiral pertur- 
bation theory we refer the reader to ||. 

4 The effective lagrangians £2? £4 and jCq 

At 0(p 2 ) the most general effective lagrangian coincides with the non-linear cr-model 
lagrangian £2 equipped with external fields ||], 

£2 = \f 2 {D^UD^ +^U + X U^) 
D^JJ = dyU - i(Vfj, + a^)U + ill (v^ - a M ), 

X = 2B{s + ip). (4.1) 

The brackets ( • • • ) denote a trace in flavour space and the unitary 2x2 matrix U 
contains the pion fields, 

U 



(4.2) 




V 



We count the Goldstone fields as O(p ), a derivative d^, the vector and axial-vector 
currents v^a^ as 0(p) and the scalar and pseudoscalar currents s,p as 0(p 2 ). The 
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lagrangian £2 contains two free parameters F and B, where F is the pion decay constant 
in the chiral limit and B is related to the order parameter ( 1 qq | ) . In terms of 
physical quantities we have 

Fir = F(l + 0(m)) , 

Ml = M 2 (l + 0(m)) , M 2 = 2Bm. (4.3) 

The transformation properties under G = SU(2)l x SU(2)r x U{\) of the several fields 
are, 

U -» g R Ug[, 
v^ + a^ -> gjtiv^ + a^Jj + ignd^ , 
v^-a^ -> - a^)^ + igLd^gl , 

s + ip -> ffi?(s + ip)^ , 

F ijjR G 5C/(2) ; = diag(0 o ,0 o ) ; G R . (4.4) 

Since the charge is not a generator of SU(2), we consider in the following the case 

(a„) =0 ; (^) ^0. (4.5) 

This condition is consistent with the transformation law ([4.4j). To calculate V^ v , we set 

3 = ml ; v M = Qw M ; p = p l r l ; a M = , (4.6) 

where Q = 1/3 diag(2, —1) is the charge matrix. 

The next-to-leading term Z4 contains two types of contributions: one-loop graphs 
generated by £2 and tree graphs involving one vertex from £4. The most general chiral 
invariant lagrangian £4 was constructed by Gasser and Leutwyler 

£4 = £ ( - 4 - ) + Cwzw , 
7 



£< 4 ) = hPi + ... (4.7) 



i=l 



where 



Pi = \{u»u,)\ _ , 
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I ^5 — ( /— /-/if 



^2 = 7 (n M n,) (uV) 

1 = i(/r^^]) . (4.8) 



The following notation was used: 



4 

lo D / 



u M = iu^DpUv) = —iuD^Wu 
X± = ut X" f ± ux^u , 



f£ = uF^v) ±u^F^u, (4.9) 
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with u 2 = U. The field strength tensors F^ U R are related to the nonabelian external fields 
r fl = v tl + a tJi ,l l j, = v f j,- a M through 



F R = d » r » ~ ®v r H 

K = d^-dJ^-i^J,} . (4.10) 

The anomaly term Cwzw contributes to V% v pi"] ; we therefore simply neglect these 
contributions. The ellipses in ( f4.7|) denote polynomials in the external fields which are 
independent of the pion variable. These terms do not contribute to S-matrix elements. 
The realization of G on the fields in Eq. (|4.9| ) is given by 



I h(4>)Ih} (4>) , (4.11) 
where the non-linear realization h(cp) defines the action of G on a coset element u{4>) via 

u(<f>) ° g R u{4>)h){4>) = h(<j>)u(<j>)gl . (4.12) 

Since the one-loop graphs generated by £2 may be ultraviolet divergent, the low-energy 
constants b L have to absorb these divergences leading to an ultraviolet finite representation 
of any contribution at 0(p 4 ). Using dimensional regularization, the low-energy constants 
are defined ||] as 

k = l r i + . i = 1, ... ,7 , 

,2u> 



A 



16tH 



{^-i(m47T + r'(l) + l)} , (4.13) 



where l\ are the coupling constants renormalized at the scale /x. The 7^ are displayed in 
Ref. §. 

Terms at 0(p 6 ) generated by Z§ contain two-loop diagrams with vertices from £2, one- 
loop graphs with one vertex from £4 and tree diagrams generated by £§. The lagrangian 
£q in the mesonic sector was recently constructed by Fearing and Scherer |}2| . In the case 
of SU(3) x SU(3), they find 111 (32) independent terms in the even (odd) intrinsic parity 
sector. In the even parity sector relevant here, the new low-energy constants that enter 
at this order have the form 

B i = B r l +Bf i = l,...,lll, (4.14) 

where the Bf contain poles and double poles at d = 4. 

The amplitude generated by tree diagrams involving one vertex from £q contributes 
only a polynomial to 77 — > ir + ir~ , which has the same structure as the divergent part in 



the two-loop amplitude. To illustrate this fact, we note that the monomials [22] 

5£ e = B 30 ( [G^] + [D^G af3 } + ) 
+ Bsi ( [G^} + [D a D^G M } + ) 

+ B 50 ( [ x ]+> ( [GH+[G>]+> (4.15) 

generate the polynomials 

8A% = [ (4B 30 + 9B 31 )s - (8B 30 + 112B 50 )M 2 } , 

«*6° = 2F^3l- (4.16) 



7 





(2) 



(3) 



Figure 1: Tree diagrams contributing to Z^. The mass in the propagators is identified with the 
physical pion mass at this order. 

One may therefore remove the divergences in the two-loop amplitude by simply dropping 
the singular part. The finite pieces one is left with may be estimated via resonance 
exchange (see below). In the following, we use the notation 



for the total contribution from Cq. 



a r { c M 2 + a% c s 

(16tt 2 F 2 ) 2 

b r ' c 
( 16vr2jF 2)2 ' 



(4.17) 



5 Results at 0(p 2 ), 0(p A ) 

The process 77 — > tt + tt~ occurs already at tree level, in contrast to the situation in the 
neutral channel, where the leading term starts at one-loop order. The amplitude at tree 
level coincides with the one obtained from scalar electrodynamics, 



A c 
B° 



1 



1 

2s 



M 2 - t 
1 



+ 



1 



+ 



M 2 — u 
1 



M 2 — t M 2 



+ 0(1), 

+ 0(1). 



(5.1) 



The pion mass in the propagators is identified with the physical one at this order in the 
expansion. The corresponding Feynman diagrams are depicted in Fig. |I[ The next-to- 
leading term Z4 contains two types of contributions: one-loop graphs generated by £2 
and tree graphs involving one vertex from £4. The term proportional to the low-energy 
constant Z3 in the lagrangian £4 is quadratic in the fields. One may therefore shift the 
pion mass squared according to M 2 — > Mq = M 2 + 2M 4l I%/F 2 and from now on we denote 
the pion propagator by (Mq — p 2 )^ 1 , also in Fig. |l[ Evaluating the several diagrams to 
0(p±) (Fig. § gives § 



A c 
B c 



1 1 

+ 



1 



M 2 -t M 2 - 

1 1 

+ 



+ 



M 2 — t Ml 



F 2 
+ 0(1). 



G(s) + 



Ia 

487T 2 



+ 0(p 2 ) 



(5.2) 
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(10) 




VWW*^ VWW*^ WA/Wj^^ wwy^ 

v-s ^ T + crossed 

VWWvJk^ VWWVPk^ WWWK^ vwO^ ,,5 ' ,9) 

(11) (12)0^ (13) (14) 



Figure 2: One-loop diagrams contributing to Z4; the graphs (4, 5) generate the finite loop function 
G(s). The graphs (10 -19) amounts to the replacement M 2 —> M% in the tree result. Graphs (8, 
9) stem from the lagrangian £4. 



The loop function G(s) is discussed in appendix [A|. The physical pion mass squared M 2 
is 



Mi = M 



32tt* 



(5.3) 



The pion decay constant F is identified at this order with the physical F n . It is a spe- 
cific feature of the process 77 — > 7r + 7r~ that, after mass renormalization, the one-loop 
contribution from £2 is ultraviolet finite. £4 generates the graphs (8,9) in Fig. |2[ This 
contribution is proportional to the finite, scale independent combination 



l A = 96TT 2 (2l r 5 -l r 6 ) =2.7. 



(5.4) 



As is discussed in more detail in section |8.4| , there is good agreement with the available 
experimental data and the chiral representation (^^). We now turn to the evaluation of 
the contributions at order p 6 . 



6 Evaluation of diagrams at 0(p 6 ) 

Up to 0(p 4 ) one has to deal with only a few Feynman diagrams for the process 77 
and the evaluation of the amplitude is straightforward. At 0(p e ) we have found remarkably 
many diagrams. To generate all graphs and to check the combinatorial factor for each of 
them, we use a beta-version of the Mathematica package FeynArts 2.0 developed by S. 



Kueblbeck and H. Eck [23], where we have included the corresponding Feynman rules 
for CHPT, adapted to 77 — > it + it~. In Fig. ^ we have plotted the full set of two-loop 
diagrams generated by £2- We have not distinguished the different intermediate pions 
running in the loops (n , it 2 , 7r 3 ). Each displayed diagram is the generating one for a class 
of subdiagrams with the same topology. 



6.1 Two-loop diagrams from £ 2 

The genuine two-loop diagrams 22 (vertex), 36,37 (box) and 45 (acnode) in Fig. || cannot 
be represented as products of one-loop integrals. A method to perform the relevant 
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A\ wwyvf^" vm/waffi wvwr^ vww*^ 



(J VWWVIL^ WWWlk^ wwwft^ VVWW*^ 

?(46) (47) Cj" (4a _50) (51 _53) (54-56) 



+ crossec 
(57-65) 



(51 -53) (54-56) 



Vfyy^ WWr^ wA-"''''' 

vwwvi^ "vODk^ wwwi^ 

(66) (67) (68) (69) ^ 

m/Jk w\rt(jQ wwvi. wvwQ) 

(70) (71) (72) (73) 

wwjw^ w^Wf^ MWflGr*^ wwO^"''' 

vwwvjt^ wwwC^ Wift/wi^ 



+ crossed 
(74-81 ) 



(85) 



W^V^ WMrtr^ WWi*^ WVW*^ 

vwvwJl w<x^ w/<x^ ^w^S^ 

/Qeh * — ' /cr?h (ooi /ont 



- crossed 
(90-97) 



Figure 3: Full set of two-loop diagrams generated by £2- The vertex (22), box (36,37) and 
acnode (45) graphs are genuine two-loop diagrams, which do not have a representation in terms 
of one-loop functions. The dash-dotted boxes indicate subdiagrams which we have treated in the 
dispersive manner. We use symmetry arguments which relate the contributions from the reducible 
diagrams (47-97) to mass- and wavefunction corrections at two-loop order. The hatched circles 
summarize selfenergy contributions to the pion propagator (graphs (b - d) in Fig. II). 
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o 2 -i 




Figure 4: Construction scheme for a subclass of two-loop diagrams. The four-point function on 
the right-hand side is the d-dimensional elastic 7T7r-amplitude at one-loop accuracy with two pions 
off-shell. The particles running between left and right are ir 1 and ir 2 . The symbol d d l stands for 
integration over internal lines with weight (6.1). 



momentum integrals is given in [||,^4,25]. The graphs (20-26, 34-41) may be generated 
according to Fig. ||, where the four-point function on the right-hand side is the d— 
dimensional elastic 7T7r-scattering amplitude at one-loop accuracy, with two pions off-shell. 
The loop integration over internal momenta with weight 

' (6-1) 



[M 2 -{l + qi ) 2 } 



is indicated by the symbol d d l in the figure. M 2 denotes the physical pion mass in one-loop 
approximation (Eq. 53) and the momenta of the pions running in the loop are (/ + qi)^ 



and (q2 — This procedure has the advantage that the nonlocal singularities generated 
by divergent subdiagrams (enclosed by dash-dotted lines in the graphs (22,36,37)) are 
automatically removed ||. 

A further comment concerns diagram (23) and (24). The tadpole in the loop propagator 
amounts to the replacement M 2 — ► M 2 if one also adds the one-loop diagram (4) in 
Fig. ^ multiplied with the wave function renormalization constant Z n given at one-loop 
accuracy. A similar argument holds in the case of diagrams (38,39,42-44). Summing up 
these diagrams allows one to introduce the physical mass in the loop propagators. 

The graphs (27-33) may be generated by multiplying the seagull diagram (Fig. |], graph 
3) with that part of the wave function renormalization constant Z^, which is generated 
by C.2- Furthermore, we use isospin symmetry arguments to handle the diagrams (47-97) 
that also incorporate the pion propagator at two-loop accuracy [24]. 

6.2 Pion propagator at 0(p 6 ) 
We compute the correlator 

5 ab ^A' c (p 2 ) = J d D xe ipx ( O|T0 a (x)/(O)|O ) (6.2) 

in the a— model parametrization Eq. ( |4.2j ). The complete propagator can be written in 
terms of the bare propagatorA'(p 2 ) = j M -i^ p i_ ie and the self energy function £(p 2 ) in the 
standard manner, 

~K(P ) = -JT2 - 2 - + -TH 2 ~ % ^P )~TH 2 - + ■■■ 

i iM^—p^—ie i M z — p z — ic i M L — p z — le 
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(a) (b) (c) (d) (e) (f) (g) 

Figure 5: Selfenergy- graphs up to 0(p e ). The contributions from diagram (a) start at 0(p 4 ), 
whereas diagrams (b-g) are two-loop corrections. The sunset diagram (c) contributes to R(p 2 ). 
The boxes in the graphs (e) and (f) denote vertices from the lagrangian £4, where the numbers 
indicate the corresponding low-energy constant. Diagram (g) contains the contribution from the 
lagrangian only. 



1 1 

i M 2 — p 2 — ie — T,(p 2 ) 

' Z +R{p 2 ) 



M 2 — p 2 — ie 



(6.3) 



In the last step we have developed the selfenergy function T,(p 2 ) in a Taylor series around 
the physical mass M 2 , 

S(p 2 ) = S(M|) - (Ml - p 2 )t(Ml) + (Ml - p 2 ) 2 S (p 2 ) , 

R(P 2 ) = = . , (6.4) 

(l + E)[l + £-(M2-p2) g (p2)] 

where Z n = (1 + S(M 2 )) -1 is the residuum of the propagator. The physical mass is defined 
as the position of the pole of the propagator and we therefore have the implicit equation 

M 2 - £(Af£) = Ml , (6.5) 

which defines the physical mass M^ order by order. 

Evaluating diagram (a) in Fig. || with the mass Mq = M 2 +2M 4: l^/F 4: in the propagator 
and solving Eq. fl6.5|) one obtains the physical pion mass up to 0(p 4 ) given in Eq. (|5.3|) . 
The wave function renormalization constant Z n and the finite part of the propagator R(s) 
become at this order 

Z 7r = l-j^F 1 (M 2 ) , R(s) = , [0(p 4 )] (6.6) 



where F\(M 2 ) is a one-loop function denned in Eq. ( |A.4| ). At 0(p®) the graphs (b-g) 
in Fig. |5| contribute to S. The diagrams (b-d) stem from the lagrangian C2 whereas 
(e) and (f) are contributions from £4. The lagrangian Cq contributes via diagram (g). 
Only the sunset graph (c) is nontrivial. The overlapping loop momenta in one of the 
propagators make the algebraic part tedious. Nevertheless for s = Ml the result of the 
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sunset selfenergy contribution may be obtained analytically. We find for the pion mass up 
to 0(p 6 ) 



Mi 



M 2 < 1 + 



M 2 / Fi(M 2 ) 



F 2 \ 2M 2 

M 4 



+ 2l 3 



F 2 {M 2 ) (f 2 (M 2 ) P m + 1 ^m] + B M 



i=l 



+ C(M 8 ) . (6.7) 



Bm collects the contributions from the Cq lagrangian and the polynomial p l M read 



PM 



PM 



Pm 



PM 



17 79 421 2 

w H u? + C(lj 3 ) 

8 12 32 V ; 

31 

14 - 15w + yw 2 + 0(u; 3 ) , 
8- 10w + llw 2 + C(w 3 ) , 
3 - 4a; + 4w 2 + C(w 3 ) . 



(6.8) 



-Fi, F2 are one-loop functions defined in Eq. ( A.4j) . For the wave function renormalization 
constant we get at 0(p 6 ) 



Z, = l~ + |f 2 (M 2 ) (f 2 (M 2 ) Pz + £ ^ + ^4 + 0(M 6 ) 



F 2 



(6.9) 



i=l 



where £?z collects the contributions from the Cq lagrangian and where 



Pz 



P l z 



P 2 z 



n 3 

Pz 



29 395 „ o N 

1 w + w 2 + C(u; 3 ) , 

6 32 v ; ' 

31 

14 - Iboo + yw 2 + C(w 3 ) , 
8 - 10a; + llw 2 + C(w 3 ) , 
2 + 0(cu 3 ) . 



(6.10) 



Furthermore R(s) becomes 

f 2 (m 2 : 



R(s) 
PR 



PR + 



1 



h F (s) 



6F 4 (s - M 2 ) 2 ' 



with 



h F (s) 
JC(x, s) 
z 2 (x,s) 



F A 

-\M 2 + (\M 2 - i,) w - (§M 2 + ^s)u 2 + 0{^) , 

rQO r 1 

/ da(3{a) / dx{ 6sM A + (1 - 12x + 18x 2 )s 2 M 2 }/C(x, s) , 

JAMi JO 



(6.11) 



1 



z 2 (x,s) (s-M 2 )x(l-x) 

m ; TT^T + 



(16^ 2 ) 2 [ z 2 (x,M 2 ) z 2 (x,M 2 ) 
M 2 - s(l - x) ] x + er(l - x) , 



(6.12) 



and where /?(c) = ^/l — AM' 2 /a. The loop function Zip stems from the sunset diagram (c) 
in Fig. || 
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(130) (131) (132) (133) (134) (135) 



Figure 6: One-loop order contributions generated by the lagrangians £2+^4 with one vertex from 
£4. The boxes denote £4- couplings, whereas the numbers indicate the corresponding low-energy 
constant. 

6.3 Contributions with one vertex from £4 

We discuss some aspects of one-loop graphs generated by £2 + £4, with one vertex from 
the lagrangian £4 depicted in Fig. ^. Only diagrams proportional to {h,hih,h) con- 
tribute to the process 77 — > iT + ir~ . The low energy constants I3 and I4 enter the two-loop 
amplitude through mass- and pion decay constant corrections. The diagrams (98,99) and 
(114,115) in Fig. ^ may be generated according to Fig. |] where the /^-vertex is contained 
in the <i-dimensional elastic 7T7r-scattering amplitude at one-loop accuracy. These graphs 
remove the sub divergences indicated in Fig. ||] (enclosed by a dash-dotted line in graphs 
(22,36) ). 

The diagrams (100,101) and (116,117) may be included in the replacement Z n seagull 
and the graphs (104-113), (120-129) are included in the replacement M 2 — > M 2 in the 
reducible Born diagrams. The remaining graphs require straightforward one-loop calcu- 
lation. 

6.4 Counterterms with one vertex from £q 

The diagrams depicted in Fig. fj] and generated by the lagrangian £§ contribute a polyno- 
mial piece only. The divergent parts of the couplings cancel the pole structure generated by 
the two-loop diagrams. We have estimated the finite parts of these couplings by resonance 
exchange (see below). 
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(136) 



(137) 



(139) 



(140) 



(141) 



(138) 

+ crossed 

WMMJ^^ VWVWPk^ (144-148) 



(142) 



(143) 



Figure 7: Tree order counterterms generated by £2 + £4 + £e with one vertex from £ 6 denoted 
by a white circle. Diagrams (137-139,141,143) denote mass corrections in the propagators. 



This concludes the discussion of the diagrams which occur at two-loop order in 77 
7r + 7r~. For a complete mathematical treatment of all genuine graphs occurring in 77 
7r+7r~(77 — > 7r°7T ) see pl,p). 



7 Checks 

The tensorial decomposition of the amplitude (2.4) leads to formfactors Cj(s, v) (see ||), 
which are linearly related through two Ward identities and through Bose symmetry, 

Li(cj) = 2c + sc 4 - (t - u)c 5 = , 

L 2 (Ci) = SC 7 - (t - U)Cg = , 

L 3 {ci) = c 5 + C7 = 0. (7.1) 

The amplitudes A° and B c may be obtained from C4 and eg, 

A = -c 4 , B = c 9 /2s. (7.2) 

We have evaluated the five formfactors Co, C4, C5, cj and eg in d dimensions and afterward 
the following consistency checks were made: 



(i) We have checked that the relations (|7.l| ) below threshold are satisfied numerically 
within machine accuracy in low dimensions. 

(ii) By letting u> — > and omitting the tree contributions generated by £2 + £4 + £-6 hi 

Fig. [?], the formfactors Cj at order have the structure 

Ci(s,u) = ^^l[F 2 {M 2 )[^ + ^ + a^ 2 + 0{^)] 

+ E lj[bl J +b\ j u; + b 2 1 u; 2 + 0(u; 3 )}} 

j=l,2,5,6 
f p(2) p(l) ^1 
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Table 1: Coefficients of the amplitude expansion in Eq. (7.3) 





i = 


i = 4 


i = 5 


i = 7 


% = 9 


a\ 


337 M 2 127 2 53 2 
216 ivl b 864* 864^ 


§M 2 - Is 


_J3_ 
432^ 


¥ 

J3_ 
432^ 


i' 

53 . 
432* 


b$ 
bf 


fM 2 s + \s 2 - \v 2 
-4M 2 s 


-f M 2 - | S 
8M 2 


-i" 




¥ 




i- 




iA2 
°0 


-3M 2 s + is 2 - \v 2 


6M 2 - is 


-I, 


i, 


is 


bf 


2M 2 s 


-4M 2 











bf 


8M 2 s + 2s 2 


-16M 2 -4s 











bf 


-8M 2 s 


16M 2 











bf 


-\M 2 s - s 2 


8M 2 + 2s 











bf 


AM 2 s 


-8M 2 












where al and 6>f , fc = 0, 1 are polynomials in the external momenta and the pion 
mass (Table |) and F 2 {M 2 ) is a loop function defined in Eq. ( |A.4| ). The coefficients 



a\ and 6 2 contain a polynomial part as well as complicated analytic functions and 
contribute to the finite part Ri of the amplitude. 

Ik) 

We have verified that the residues P£ are polynomials in the external momenta 



and the pion mass. These polynomials obey Eq. (7.1) as well, 

L m (if°) = 0, m= 1,2,3, A; = 1,2. (7.4) 

We have checked numerically that 

L m (Ri)=0, m = 1,2,3, (7.5) 

within machine accuracy below threshold. To go to the physical region for pion- 
pair production we have written fixed-i dispersion relations for the box- and vertex 
parts. Matching the dispersive analysis with the non-dispersive expressions below 
threshold ensures that L m {R,i) = is also satisfied above threshold. 

(iii) We have worked out the S-wave projection h% + (s) ( ^S| ) for the helicity amplitude 
H9, and we have verified numerically that these amplitudes admit the correct phase 
at s > 4M 2 , given by the (tree+one-loop) elastic tttt scattering S-wave phase shifts 
(in the appropriate isospin decomposition). 



8 Amplitude and cross section to two loops 
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8.1 Analytic results of the amplitude 

We find for the two-loop amplitude A° 

A c = A c 2 +A ( i + AC + 0(p i ) , (8.1) 

or 

A ° - { nh + + h [rMs) + ^' + u ° + Pi + ° (pi) ■ <8 - 2) 

The unitary part U% contains s, t and u- channel cuts, and P% is a linear polynomial in 
s. Explicitly we find, 

U a = ^G(s) [ (2M% - 4M*a + 3s 2 )J(s) + C c (s, k) ) + ^^rM 8 ) 

+ 288^4 ( * " AM -^ S ^ + 4 [ s0{s) + 2M *^ (s) " 3 =J (s)) ] d ™ } 

+ ( s ~ ± M l){H{s) + 4 [ sG{s) + 2M 2 (G (s) - 3 J (s)) ] d 2 00 } 

+ A c A (s,t,u), (8.3) 

with 

-12M% + YIMlsh ~ 12M 2 s + 12M 4 } , 
d 2 = I(3cos0 2 -1). (8.4) 

The loop functions J etc. are displayed in appendix [A] and G n (s) in Eq. ( |8.2| ) stands for 
G(s) evaluated with the physical mass. The term proportional to (igg in U° contributes 
to D-waves only. For A A see below. The polynomial part is 

P a = (16 ^ F 2)2 WiMl+als] , 

1 337 5 

°i = a i r + tt{ 4/2 + *(-10ii + 18Z 2 - 12f A + — ) - -h - 5f a + I2/4/A + 4} , 
y bo 

„ 1 , , X — 3— „ — 127, 5— 5— „ — 21-, 

«2 = a c 2 ' r --{/ 2 +/(-Z 1 + -/ 2 + 3/A + ^ r )-^Zi- i Z 2 + 3/ A + y}, 

M 2 

I = ln^f. (8.5) 
The result for B c reads 

B c = BC + BC + 0(p 2 ), (8.6) 

or 

with the unitary part 

U C B = ^^pi i\ih ~\) + (h ~ l)}H(s) + A c B (s,t,u) . (8.8) 
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For the polynomial part we find 



(16^2)2 ' 

b c,r _ J_n2 
18 1 



53, 
24' 



12 



I* 2 + l } ■ 



(8.9) 



The integrals B (s,t,u) contain contributions from the two-loop box-, vertex- and 
acnode graphs and also from the reducible diagrams manifest in the function R(s) given 



in Eq. ( 6.11 ). We discuss the size of these corrections in section 8.3 . 



8.2 Low— energy constants 

The Born contribution to the charged pion-pair production contains M n as the only 
parameter. At next-to-leading order three parameters occur, namely, M^F^ls — h- 
Going to the two-loop correction the number of parameters increases. Once the program 
described in section ^ is carried through, the ultraviolet finite and scale independent 
amplitudes A c and B c contain the set M n , F;l\ ,l r a (i = 1, 2, 3, 5, 6); a{ c , dj and b r,c of 
low-energy constants. The couplings stem from the O{oj)— part in the expansion of the 
£4 low-energy parameters, whereas d{ c a 2 ' c and b r,c are couplings from the lagrangian £,%. 
F is related to the pion decay constant F n M 



F W = F 



1 + 



M2 

F 2 



1 



16vr 2 



In 



M2 



+ 0(M 4 ) 



.10) 



We may therefore replace F by F n at the expense of introducing l\. The expressions for 
the loop-amplitudes simplify if one uses the scale-independent parameters k instead of l\ 
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32tt 2 



(li + In — 5- 
i 1 



(8.11) 



The values of the \ are depicted in Table ^ together with the experimental sources and the 
values for the 7$. The couplings contribute a polynomial piece only and may therefore 
be absorbed into the low-energy constants at order p e [27|. We are therefore left with 



1 ! 



a 2 ' and b r ' c as the only new unknowns. The coupling constants occurring at 0(p ) in 
the low-energy expansion are dominated by the low-lying vector, axial-vector, scalar and 
pseudoscalar resonances in such a way, that they practically saturate the corresponding 
couplings of the 0(p 4 ) effective lagrangian j|, 28 1. 

Here we assume that this is also true at order p 6 . We estimate the renormalized 
couplings a^' c ,02' c and b r ' c by replacing them at a scale \i = 500MeV . . . lGeV by the 
contribution from resonance exchange. Let 



£ / R + / r Gu), I = a h a 2 ,b, 



$.12) 



R=p,a\,b\ 



V 



(J PC = 1 ) and axial-vector ( J PC = 
) exchange [11, 2S-32| with meson mass Mr < 1.2 GeV. The non-resonance 



where the sum denotes contributions from vector 
1++ 



contribution of these low-energy constants has to be chosen in such a way that the final 
result of the amplitudes A c and B c is scale independent. Our estimate for I c ' r (M p ) 
consists in setting I r (M p ) = 0. 
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Table 2: Low-energy constants and sources. In the fourth column we have plotted the % evaluated 
in gf§ 



i 


k 


Source 


7i 


1 


-1.7 ± 1.0 


K e 4 , 7T7T — > 7T7T 


1/3 


2 


6.1 ±0.5 


K e £ , 7T7T — > 7T7T 


2/3 


3 


2.9 ±2.4 


S 1 ?/ (3) mass formulae 


-1/2 


4 


4.3 ±0.9 


Fk/Fk 


2 


5 


13.8 ± 1.3 


it — > ei^7 


-1/6 


6 


16.5 ± 1.1 


< r 2 >^ 


-1/3 



Table 3: Resonance contributions to the constants a r { c , a^ and b r,c . The calculation is done 
in [ p4[ , the uncertainties are more generous than quoted in that reference. 



jr 




E R i R 


P 


ai 


h 


r,c 


-3.28 








-3.3 ±1.65 


r,c 

a 2 


1.23 


-0.35 


-0.13 


0.75 ±0.65 


b r ' c 


0.20 


0.18 


0.06 


0.45 ±0.15 



The quantities I R are evaluated in |^4|. The results of this calculation are displayed 
in Table || where the individual resonance contributions I R are listed. Column 5 contains 
the sums of the several contributions provided with a systematic error. We associate a 
50% uncertainty to the contributions generated by (axial-) vector exchange. 

It is useful to consider furthermore the helicity amplitudes and H+_ and the 

corresponding low-energy constants /i± c and h r s ,c , 

TjC,2loops _ ti+ c M 2 + h r s ,c s 
++ ~ (16vr 2 F 2 ) 2 + "' ' 

rrC,2loops _ 8(M 4 — tu) rjC 
+~ - s ( 167r 2 F 2)2 n - ' 

hl c = a r { c + 8b r ' c , h r f = a r i c - 2b r > c , h r f = b r > c . (8.13) 

Adding the errors quoted in Table |3| in quadrature, we find 

h r _f{M p ) = 0.3 ± 2.0 , 
h r s ' c (M p ) = -0.15 ±0.7, 

h r l c {M p ) = 0.45 ±0.15. (8.14) 
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This completes the determination of the parameters which occur at two-loop order in 

77 — ► 7T + 7T~. 



8.3 Numeric results of the amplitude 

The characteristics of the two-loop corrections may be seen in Fig. H where we have 
plotted the real and imaginary part of the helicity amplitudes W^^M^H^ ( |2.10 ) at 
t = u. All curves start at E = 280 MeV and the crosses refer to the center-of-mass energy 
of the 7r + 7r~ system in 100 MeV steps []. The solid line incorporates all contributions 
stemming from the one- and two-loop calculation only. The dash-dotted line is the 
same amplitude without A^ B . The dotted line corresponds to the one-loop result (apart 
from mass renormalization terms, H9__ does not receive any contribution from one-loop 
diagrams). To get some feeling about the size of the genuine two-loop diagrams, we have 
drawn with a dashed line the contributions from the sum of one-loop and acnode graphs, 
with a dashed-double dotted line the sum of the one-loop and box diagrams at the scale 
H = 770 MeV in the MS chiral scheme. 



8.4 Cross section 77 — ► n + n 

The cross section for 77 — ► tt + tt~ in the framework of SU{3) x SU(3) CHPT was worked 
out up to 0(p 4 ) in ||. The amplitude at this order contains only the three low-energy 
parameters F n , M n and Z a • We note that the increase in the cross section around the 
peak at = 306MeV due to the one-loop correction is ~ 15% where the renormaliza- 
tion scale independent counterterm contribution ~ l/\ amounts ~ 50% of this correction. 
Furthermore it is instructive to compare the cross section with the one from 77 — > 7r°7r° 
scattering ||. In the neutral case the size of the cross section is two orders of magnitude 
smaller than in the charged channel, since the neutral amplitude starts out with one-loop 
diagrams, whereas in the charged case the Born amplitude contributes dominantly to the 
cross section: At low energies the photon couples to the charge of the pions. As a result 
of this, the 77 — ► 7r 7r° cross section is small, while that for 77 — ► -k + it~ is large. 

The plot in Fig. ^ shows the total cross section a(s ; | cos 9 \ < Z = 0.6) as a function 
of the center-of-mass energy E = yfs. Here we choose Z = 0.6 to compare with the 
available data from the Mark II collaboration The Born contribution is displayed 

with a dotted line whereas the born+one-loop cross section is shown with a dashed line. 
Also at this order the chiral expansion is in remarkable agreement with the available data. 
The solid line incorporates all contribution up to 0(p G ). It is seen that the two-loop effect 
changes the one-loop result only very little. Switching off the integrals B does not 
affect the cross section visibly. This fact is also clear by consulting Fig. ||, from where 
it is seen that our explicit, compact representation of the two-loop amplitudes is rather 
accurate. We note that the two-loop correction, mainly dominated by the effect of final 
state 7T7T interaction, enhances the cross section very close to threshold and then suppress 
it above E = 350 MeV, see Fig. |9| Finally, the dashed-double dotted line displays the 
result of a dispersive analysis (Fig. 7 in Ref. p| ). In that calculation, use was made of a 
doubly subtracted dispersion relation with the subtraction constants determined in terms 
of chiral counterterms, incorporating also heavy meson exchange. 

3 We use F w = 92.4 MeV, M n = 139.6 MeV and h = -1.7 ± 1.0, h = 6.1 ± 0.5, l 3 = 2.9 ± 2.4, 
h = 4.3 ± 0.9, Ta = 2.7 ± 0.4 throughout. 
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Figure 8: Real and imaginary part of the helicity amplitude 10 -3 M%H+± at t = u. The solid line 
incorporates all contributions from one- and two-loop corrections. The dash-dotted line is the 
same amplitude with the omission of B . The dotted line corresponds to the one-loop result 
(apart from mass renormalization terms, H+_ docs not receive any contribution from one-loop 
diagrams - the dotted line is therefore absent in the lower figure). The dashed line is the sum of 
the one-loop and acnode, the dashed-double dotted the sum of the one-loop and box contribution 
at the scale /i = 770 MeV in the MS chiral scheme. Finally the crosses refer to the center-of-mass 
energy of the tt + tt~ system in 100 MeV steps. 
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Figure 9: The 77 — > ir + ir~ cross section er(s ; | cos 6* | < Z = 0.6) as a function of the center-of- 
mass energy E, together with the data from the Mark II collaboration JlOj ] . We have added in 
quadrature the tabulated statistical and systematical errors. In addition, there is an overall nor- 
malization uncertainty of 7% in the data |l0[ . The solid line is the full two-loop result, the dashed 
line corresponds to the one-loop approximation j^] and the dotted line is the Born contribution. 
The dashed-double dotted line is the result of a dispersive calculation performed by Donoghue and 
Holstein (Fig. 7 in Ref. Q). 



The two-loop result thus agrees very well with the data and shows furthermore a very 
good agreement for ^fs > 350 MeV with the dispersive analysis of Donoghue and Holstein. 

In Fig. 10 the effect of the low-energy constants is shown. The solid line corresponds 
to the two-loop result evaluated at B = and without resonance exchange. (It turns 
out that the contributions from the resonances to the cross section are negligible below 
450 MeV.) The dashed line corresponds to li = 0, and the dash-dotted line is obtained by 
setting l\ = I3 = 0. Therefore the increase in the cross section is due to l2,h and Ia, as 
was also shown in the neutral case Hi. 



9 Compton scattering and pion polarizability 

Compton scattering on pions jtt + —* jtt + is related to piompair production by crossing 
(Pi;<?i) At — * ~ (Pi;9i) At ; s ^ The corresponding scattering amplitudes may again be 
expressed in terms of the amplitudes A c and B c as functions of the Mandelstam variables 
in the Compton channel, 

* = (Qi +Pi? , t = (Q2 ~ qif ,u = (p 2 - qi f . (9.1) 

In the present section we discuss shortly the Compton cross section. Afterwards we work 
out the chiral expansion of the charged pion polarizabilities at next-to-leading order. 
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Figure 10: The 77 — > tt + tt cross section a(s ; | cos 8 | < Z = 0.6) as a function of the center-of- 
mass energy E. The solid line is the two-loop result evaluated at B = and without resonance 
exchange. The dashed line corresponds to k = 0, whereas the dash-dotted line is obtained by 
setting £1 = £3 = 0. 

9.1 Compton cross section jn + — > 77r + 



The total cross section at tree level is given by 33] 

Ml + s 



' (s) = 4a 2 ir 



s(M2 - s) 3 



M^-s 2 + 2Mls\n^ i 



(9.2) 



This quantity is plotted in Fig. ^ as a function of the center-of-mass energy E 17T . At 
threshold one obtains the total Thomson cross section 

a^ + ^\s = Ml) = ^rl, (9.3) 

where r e = a/M^ is the classical electromagnetic radius of the pion. Comparing the 
Compton cross section with the pion production case we note that the former is one order 
of magnitude bigger. The increase in the cross section due to one-loop corrections is 
smaller than one percent and the two-loop result is again smaller by a factor of 3 . As a 
result of this, corrections in the charged Compton channel are negligible. 

In the neutral channel, the leading term is generated by one-loop graphs. The cross 
section at E = 350 MeV is three orders of magnitude smaller than in the charged case 
and tends to zero at threshold. Two-loop corrections are substantial, and the result for 
the Compton amplitude is therefore not very reliable at this order |J. 
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Figure 11: The Compton cross section as a function of the center-of-mass energy E = in 
Born approximation. One- and two-loop corrections are tiny in this case and would not be visible 
in the figure. 



9.2 Chiral expansion of the pion polarizabilities 

Among other fundamental parameters (mass, charge, magnetic moment, ...) the elec- 
tromagnetic polarizabilities characterize a composite system like a hadron |34]]. They 
parametrize the first correction in the Taylor series expansion of the Compton amplitude 
in photon energies at threshold, 



T 



ei • e 2 *(e 2 - inM^LO^) - 4vrM^(gi x ?i) • (q 2 x e 2 *) + . . . , (9.4) 



with q 1 ^ = (u>i, qi). In terms of the helicity amplitudes H+ + (J = 0) and H+_(J = 2) one 
has 

— OL — 

a~± ± B~± = — — lim lim , (9.5) 

indicating that the combination a n ± — f3 n ± is pure S-wave while a n ± + f3 n ± is pure D-wave 
in the 77r ± — > 'yir^ channel. Below we also use the notation 

(q ± f3) c = a n ± ± p n ± , 

(a ± (3) N = a n o ± P w o . (9.6) 



From the Compton scattering amplitude at threshold Eq. (9.4) expressed through the 
functions A° and B° (see Eqs. Q8.2| , |8.7| )) it is straightforward to determine the chiral 
expansion of the pion polarizabilities in terms of the quark masses to two loops, 

«* ± «"' c = i^kn Ar + w% + ° (M ' 4) ' (9J) 
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The leading result is given in jlj| 

•4 I I -V± 



"3 / \ 3 





2 Ja 



(9.8) 



The next-to-leading order terms B± have been determined in ||. We find the following 
expressions for the coefficients B± 



a? = ^CM)-|{i(i + ir 1 + |&)-H I+ ir 1 + |& + § + A 

M 2 

Z = ln^f, (9.9) 



with 



A'; = -18(16tt 2 F 2 ) 2 lim lim AS = -8.69 , 

3 



, r9 (167TW lim lim (AS + 8Mf A§) = -8.73 . (9.10) 

Our final result of the charged pion polarizabilities including corrections up to 0(p e ) reads 
(ji = M p ) 

{a + (3) c = 0.3 ±0.1 (0.0), 
(a-0) c = 4.4 ±1.0 (5.4 ±0.8), 
a w ± = 2.4 ±0.5 (2.7 ±0.4), 

f3 w ± = -2.1 ±0.5 (-2.7 ±0.4). (9.11) 

The numbers in brackets denote the leading order result. The estimate of the errors 
only stems from the uncertainties in the couplings and do not contain effects from higher 
orders in the quark mass expansion or any correlations. (The errors are more generous 
than quoted in Ref. p4l].) We use furthermore l/\ = 2.7 ± 0.4 Q extracted from the tree 
result (0(p 4 )) of the radiative pion beta decay. Note however, that corrections in this 
process may change the error accordingly. 

9.3 Data on pion polarizabilities 

Performing Compton scattering on pions is not easy. Nevertheless, it is possible to probe 



the pion polarizability according to Fig. 12 by measuring the Compton scattering am- 



plitude: (a) Pion photoproduction in photon-nucleus scattering 7p — ► 77r + n ]17j, (b) 
Radiative pion nucleus scattering (Primakoff effect) ir~ — Z — » 7r _ — 7Z [^] or (c,d) Pion- 
pair production in e + e~- collisions. Analyzing the data with the constraint (a + (3) c = 
gives 

' 40 ± 24 (Lebedev [171 ) , 
(a-(3f = { (9.12) 
13.6 ± 2.8 (Serpukhov |I| ) . 
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Figure 12: Compton scattering 777 — ► ^tt from the charged pion has been performed both via 
pion-photon production (a) and via radiative pion scattering (b) from a heavy target. Figures (c) 
and (d) show photon annihilation 77 — > irn processes. The double lines denote nucleons. 



Relaxing the constraint (a + (3)° = 0, the Serpukhov data yield |55] 



(a + /3) c 
{a-(3f 



1.4 ±3.1 (stat.) ±2.5 (sys.) , 
15.6 ±6.4 (stat.) ±4.4 (sys.) , 



(9.13) 



where we have evaluated (a — (3) c from (3^ and (a + f3) c as given in Ref. [35|, adding the 
errors in quadrature. 

The amplitude 77 — ► 7r + 7r~ at low energies is mainly sensitive to S-wave scattering. 



In Ref. [36| unitarized S-wave amplitudes have been constructed which contain (a — (3) c 
as an adjustable parameter. A fit to Mark II data [pi] gives 



(a - PY 



4.8 ± 1.0 



(9.14) 



The result ( |9.14 ) contradicts the Serpukhov analysis Eq. ( p. 12 ), Taking into account also 
D-waves, Kaloshin et al. [37] find 



0.22 ± 0.06 ( Mark II [10] ) 
0.30 ± 0.04 ( CELLO 



(9.15) 



A detailed analysis of the same data has also been performed in |39j. The authors conclude 
that the errors quoted in Eq. ( p,15| ) are underestimated, see also Ref. [40|. By reanalyzing 
the present data on the angular distribution of the 77 — > tt + tt^ reaction with using of the 



unitary model for helicity 2 amplitudes, Kaloshin et al. [41] find 



0.22 ± 0.07 (stat.) ± 0.04 (sys.) ( Mark II [10] 
0.33 ± 0.06 (stat.) ± 0.01 (sys.) ( CELLO p 



(9.16) 



(We have taken into account that the definition of the polarizability in Refs. [36,37,41] is 
47T larger than the one used here). 

Comparing with the data, the two-loop result for (a + 0)° agrees within the error 



bars with the one found by Kaloshin et al. [[37],|4l[] (note the remark after Eq. ( 9.15 )) 
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and also with the data from Serpukhov |35|]. The value for (a — (5)° agrees within 1^ 
standard deviation with the measurement done at Lebedev [JT^] and also with the data 
analysis using unitarized S-waves in [36], however, appears to be inconsistent with the 
value quoted in fD| . 

For a more detailed discussion of the issues considered in this section, we refer the 
reader to Ref. [42|. In the following we shortly comment on the relation between the 
polarizabilities and the crossed channel reaction 77 — ► 7r + 7r~. 



9.4 Pion polarizabilities from 77 — > 7i + n data 

One may attempt to extract the polarizabilities from 77 — > ir + ir~ data in the following 
manner (T^j. Our chiral representation at 0(p 6 ) contains three new low-energy constants 
/ij! c , h r s > c . The two parameters \i± may be replaced by the polarizabilities (q=f/?) c , whereas 
h r s ' c may be determined e.g. from resonance exchange. In Fig. ^ we have plotted the 
77 — > ir + ir~ cross section for a fixed value (a + [3) c = 0.31 and a fixed value h r s ' c = —0.15, 
varying (a — 0) c between 2.21 and 8.84. The sensitivity of the cross section to a change 
in (a + (3) c is even weaker. Even a 100% change in (a ± (3) c is also consistent with the 
low-energy data. See also Fig. 9 in Ref. [11]. 



On the other hand by choosing a value for E and a one may extract for a given value 



h r s ,c the relation between (a — (3) c and (a + (3)° . This is plotted in Fig. 14 for different 
values of hg' e , where one can see that the slope of (a — (3) c is small and also changes in 
h r s ' c do not affect (a — 0)° significantly. We note that the sensivity of (a — f3) c due to 
variation in a or in the low-energy couplings li is very high. A one percent change e.g. in 
the cross section induces a 12% change in (a — j3) c ' . 

We confirm the statement [jll] that, using the chiral representation of the 77 — > 7r + 7r~ 
amplitude as an interpolation, the statistical uncertainty of the present low-energy data 
as well as the uncertainties in the low-energy couplings do not allow one to pin down the 
charged pion polarizabilities to a reasonable accuracy. The same is true in the neutral 
channel ]S, 11, 43[. 



10 Summary and conclusions 

(i) We use SU{2) x SU{2) x U{\) chiral perturbation theory to investigate the 77 — ► 
7r + 7r~ process up to two loops. We found that far more than 100 different diagrams 
contribute at C(p 6 ). Among these are genuine two-loop diagrams that cannot be 
represented as products of one-loop contributions (box, vertex, acnode and sunset). 
We find it important to take these diagrams into account, because they contribute 
substantially even near the threshold. The three new low-energy constants appearing 
at 0(p 6 ) have been estimated via resonance saturation by including vector- and 
axial-vector mesons (J PC = 1— ,1+-,1++). 

(ii) The one-loop approximation of the amplitude reproduces the available data on 
charged pion pair production from the MARK II experiment remarkably well. The 
corrections at 0(p e ) induce further small changes in the cross section, that then 
agrees rather well with the dispersion theoretic calculations j}Tl] , |l2"| at and also sub- 
stantially above the threshold region. The enhancement (suppression) of the cross 
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Figure 13: The two-loop result for the 77 — * ir + ir~ cross section parametrized through the 
polarizabilities. We set (a+(3) c — 0.31 and hg' c = —0.15, varying (a— /3) c between (a — (3) c = 2.21 
(dash-dotted line) and (a — (3)° = 8.84 (dashed line). The solid line corresponds to (a— (3) c = 4.42. 




Figure 14: Relation between (a + 0) c and (a — (3)° by fixing a point in the cr-E-plane (E=346 
MeV, cr=223.42 nb), and for different values of h r s ' c . The solid line corresponds to the central value 
h r s ' c — —0.15 whereas the upper (lower) line is evaluated with h r s ' c = —0.55(0.25). 
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section very close to threshold (above E = 350 MeV) is mainly due to tttt rescat- 
tering and renormalization of the pion decay constant, see Fig. |9[ The couplings 
\i± and h r s ' c contribute with a negligible amount below E = 450 MeV, they become, 
however, more important at higher energies. 

(iii) The cross section for jtt^ 1 — > 77r ± is fully dominated by the Born terms. One- and 
two-loop corrections are completely negligible in this case. 

(iv) The expansion of the amplitude at the Compton threshold in powers of the photon 
energies allows one to extract the pion polarizabilities in powers of the quark masses. 
We note that the sum of the polarizabilities is positive, as is suggested by the unsub- 
tracted forward dispersion relation (see e.g. p4[ ). Furthermore, (a + f3) c vanishes 
in the chiral limit m u = rrid = 0. The chiral expansion of the polarizabilities con- 
tain chiral logarithms ~ M n In 2 M n and ~ M n In M n which are suppressed by small 
Clebsch-Gordan coefficients. Nevertheless their effect on the pion polarizabilities is 
substantial. The effect of the low-energy constants h+ c on the pion polarizability is 
large (~ 50% of the two-loop result). It is therefore essential to make additional ef- 
forts in estimating the 0(p e ) low-energy constants in order to reduce the theoretical 
uncertainties in the polarizabilities. 

(v) We comment on the comparison of the chiral predictions for (a±/3) c with the data. 

- Our result (a — j3) c = 4.4 ± 1.0 includes the leading and next-to-leading order 
terms. It agrees within 1^ standard deviation with the result 40 ± 24 found at 
Lebedev fl7|. On the other hand, it is inconsistent with the value 13.6 ± 2.8 



determined at Serpukhov 18 



The analysis done by Kaloshin et al. p6, 37,41] for (a ± j3) c agrees within the 



error bars with the chiral predictions. Beware, however, the remark after Eq. 

(EHD- 

The value {a+(5) c = 0.3±0.1 includes the leading order term, generated by two- 



loop graphs. It is in good agreement with the result {a + f3) c = 0.39 ±0.04 2£], 



obtained from a forward angle dispersion sum rule based on the optical theorem 
and evaluated in a model dependent way. (Note, e.g., that the lagrangian used 
in this model is not chiral invariant.) 

(vi) In principle, the chiral representation of the 77 — ► 7r + 7r~ amplitude may be used as an 
interpolation to extract the pion polarizabilities from the low energy data. However, 
the cross section is rather insensitive to (a±/3) c , and a reliable determination seems 
not to be possible in this manner |0,E2|. 

(vii) In order to clarify the experimental situation, new experiments to determine the pion 
polarizabilities have been planned at Fermilab (E781 SELEX), Frascati (DA<£NE), 
Grenoble (Graal facility) and at Mainz (MAMI). We refer the reader to the section 
on hadron polarizabilities in Ref. [2C] for details. 
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A One— loop integrals 

We use dimensional regularization and set 

oj = | - 2 < 0. (A.l) 
Physical results will be obtained by letting d — > 4 (u> — > 0). 

A.l The integrals F m (z) 

Most of the two-loop graphs can be expressed through the one-loop function F rn (z). We 
use the following representation : 

r d d l l 

The F m satisfy the recursion relation 



and are given by 



qnri j hi 

zF m (z) = — F ro _i(z) , m>2, (A.3) 

m — \ 



F m (,) = (4vr)- 2 - r(m 2 ") ^+2-m , ro >i. (A.4) 



A. 2 The loop integrals J(s), J(s) and J (s) 

We denote by J{s) the standard one-loop integral (s = p 2 ) 

M /■ ^ _J 1 , . 

{S) J i(2n) d M 2 -PM 2 -(l+ p) 2 • 1 ] 

After Feynman-parametrization one gets 

J(s) = f 1 dxF 2 {z) , 
Jo 

z = M 2 - sx{l - x) . (A.6) 

J(s) is real analytic in the complex s-plane, cut along the positive real axis for IZes > AM 2 . 
The pole at u = may be isolated by writing 

J(s) = J(0) + J(s) , (A.7) 

where J(0) contains the pole, 

J(0) = F 2 (M 2 ) = -t^2 - + 0(1) ; w - , (A.8) 
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and J(s) is the finite loop integral 



J » = -i^/„' 

At small s, after expanding the integrand one finds 



dxln 1 — tt^x(1 — x) 
' M 2 ' 



1 00 / <j \ n 



n 2 



16vr 2 ^ \M 2 J n(2n+l)! ' 
For s > AM 2 the logarithm generates an imaginary part 



16vr 



Explicitly, 



16vr 2 J(s) = < 



JmJ(s) 

/3 (l n i_| + ivr ) + 2 
2 _ 2(^^)1 arctan ( 



s > AM' 



AM 2 



4M T '- 



0hi 



0-1 
13+1 



+ 2 



In the text we also need 



AM 2 < s 
< s < AM 2 
s < . 



J (a) = J(a) - sJ'(O) , 



(A.9) 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 



A. 3 The loop integrals G(s),G (s) and H(s) 

The integral G(s) may be obtained from the tensorial integral 



d d l 
i(2tt) c 



{2l + qi y{2l-q 2 y 



D X D 2 



<T ~ -<£<&) G{s) , (A.14) 



where 



D = M 2 -l 2 ; D 1 = M 2 - (I + q{) 2 ■ D 2 = M 2 - (I - q 2 ) 2 , 



(A.15) 



G is analytic in the complex s- plane, cut along the positive real axis for IZes > AM 2 . At 
small s, after expanding the integrand one finds 



G(s) = 



1 



E 



r - y 



(n!) 5 



l^^KM 2 ) (n + l)(2n+l)! " 



(A.16) 
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For s > AM 2 the logarithm generates an imaginary part 



8stt 1 1 - /3(«) 



JmG(s) = ^lnl±±m\ , s>AM 2 



1 



4M 2 



(A.17) 



Use of the Spence function Li 2 , especially of 

-y x 



Li 2 (y) + Li 2 { 



Li 2 (y) 



^ln\l-y), 

I" - ln(l - x) , 
Jo x 



gives 



16vr 2 G(s) 



In the text we also need 



1 + ^_ ( In f^§ + in 



AM' 2 



arctan 2 ( 



s 

4M^ 



1 + Ml j n 2 §=± 



= G(s)-sG'(0) 
sG'(O) ' " 



4M 2 < s 
< s < 4M 2 
s < . 



192vr 2 M 2 ' 

The loop function H (s) is defined in terms of G and J, 

H{s) = {s- WM 2 )J(s) + 6M 2 G(s) . 



(A.18) 



(A.19) 



(A.20) 



(A.21) 
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